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Photon-photon scattering, due to photons interacting with virtual electron-positron pairs, is an 
intriguing deviation from classical electromagnetism predicted by quantum electrodynamics (QED). 
Apart from being of fundamental interest in itself, collisions between photons are believed to be 
of importance in the vicinity of magnetars, in the present generation intense lasers, and in intense 
laser-plasma/matter interactions; the latter recreating astrophysical conditions in the laboratory. 
We show that an intense photon pulse propagating through a radiation gas can self-focus, and 
under certain circumstances collapse. This is due to the response of the radiation background, 
creating a potential well in which the pulse gets trapped, giving rise to photonic solitary structures. 
When the radiation gas intensity has reached its peak values, the gas releases part of its energy 
into ‘photon wedges’, similar to Cherenkov radiation. The results should be of importance for the 
present generation of intense lasers and for the understanding of localized gamma ray bursts in 
astrophysical environments. They could furthermore test the predictions of QED, and give means 
to create ultra-intense photonic pulses. 
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In classical electrodynamics, as described by the 
Maxwell equations, photons do not interact as long as 
there is no material medium present. However, due to 
the interaction of photons with virtual electron-positron 
pairs, quantum electrodynamics (QED) predicts photon- 
photon scattering in vacuum Q . This is commonly mod¬ 
eled by the Heisenberg-Euler (H-E) Lagrangian, which 
neglects dispersive effects. The H-E Lagrangian gives 
rise to cubic nonlinear corrections to Maxwell’s vacuum 
equations, similar to the self-interaction terms encoun¬ 
tered in optics of Kerr media |2|- The H-E corrections 
give rise to both single particle effects, such as photon 
flitting iH, lensing effects in strong magnetic helds 
[ 3 , like the ones in magnetar environments |^, and to 
coherent held effects such as harmonic generation Q and 
self-focusing of photon beams ^ . Efforts to detect these 
collisions are being made by usiM state-of-the-art super¬ 
conducting microwave facilities . Recently, it has been 
shown theoretically that QED effects can give rise to 
two-dimensional collapsing photonic structures in a ra¬ 
diation gas 0 , which could be of importance for photon 
propagation in stellar atmospheres and in the early Uni¬ 
verse. Studies of photon pulses in a radiation background 
and of optical pulses in nonlinear media reveal that they 
share common features which can be described math¬ 
ematically by a Schrodinger equation with a nonlinear 
potential [MQ- 

Dispersive effects can play an important role for short 
optical pulses when the spatial gradients and time vari¬ 
ations become large. One of the most important generic 
effects of the dispersion is to permit pulse splitting along 
the direction of propagation [iflllllj which has also been 
experimentally verified M The pulse splitting is of 
great interest in applications to normal dispersive media, 
where the collapse of light pulses can be arrested (which 
is not possible in anomalous dispersive media 


ing rise to a train of high-intensity pulses; these pulses 
may then work as a source of white light generation |l4| . 
The quantum vacuum also posses such dispersive effects 
[13, which can be of importance in the present gener¬ 
ation intense lasers |Hi| where ultra-short high-intensity 
photon pulses will be produced. Rapidly varying fields 
are also of interest for the frequency up-shift of photons 
in photon acceleration |l7l IlSl |l9j , which is an impor¬ 
tant ingredient in studies of plasma-based charged parti¬ 
cle accelerators and laser-plasma/rnatter interactions 
where ultra-strong helds |2l[l22L l2J|are expected to reach 
from peta- to zeta-watt powers j24| . In this Letter, we 
present for the hrst time results dealing with the non¬ 
linear propagation of three-dimensional intense photon 
pulses in vacuum where QED effects play a major role. 
Specihcally, we report on new features of nonlinear prop¬ 
agation of a linearly polarized intense photon pulse on a 
radiation gas background, assuming that there is no pair 
creation and that the held strength E is below the crit¬ 
ical Schwinger held, i.e. uj <C mf,c^/h ~ 8 x 10^°rads“^ 
and \E\ <C meC^/eXc — 10^®Vm“^, respectively, where 
uj is the photon frequency, uie the electron mass, e the 
magnitude of the electron charge, and Ac the Comp¬ 
ton wavelength. The derivative corrections to the H-E 
Lagrangian gives rise to a nonlinear dependence of the 
photon frequency on the wavenumber, which is shown 
to permit three-dimensional self-focusing. For moder¬ 
ate intensities of the photon pulse, the self-focusing will 
be followed by pulse splitting and the formation of sta¬ 
ble photonic solitary-pulses. However, for high initial 
powers of the photon pulse, the latter undergoes catas¬ 
trophic collapse, giving rise to field amplitudes exceed¬ 
ing the Schwinger limit 10^®Wcm“^. At these inten¬ 
sities, our theory breaks down and one must consider 
higher-order nonlinear effects and pair productions. It 
is argued that higher-order nonlinear effects become im- 
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portant before the Schwinger limit is reached, arresting 
the collapse and giving rise to ultra-high intensity three- 
dimensional solitary photonic pulses. Thus, the results 
presented here gives, apart from its immediate fundamen¬ 
tal interest and astrophysical applications, a mechanism 
for creating ultra-high intensity photon pulses. 

The evolution of an intense short photon pulse and the 
radiation background is g overned by the Karpman-like 
system of equations P ll5[ l2^E3 | 
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where <^o and (0 is the background and perturba¬ 
tion of the radiation energy density, respectively, Vg = 
c{l — fi — ^^5) is the group velocity, fdz — 2/r^(5 the 
vacuum dispersion coefficient, and k ~ cfco(l + fJ.S)jS'c 
the vacuum nonlinear refraction parameter Here, 
d = fcg/^c and /r = and the critical parameters 

k~^ ~ 10“^^ m, the Compton wavelength divided by 27r, 
and S'c ~ 10^^ Jm“^ are defined by the QED properties 
of the vacuum. 

If the time response of the radiation background is 
slow, then Eq. may be integrated to yield S’ ~ 

3/reol^'P and from Eq. llT^ we obtain the standard equa¬ 
tion for analyzing ultra-short intens e p ulses in normal 
dispersive media, see Refs. i El Ell l and references 
therein. It is well-known that the evolution of a pulse 
within this equation is modulationally unstable d, and 
displays first self-focusing, then pulse splitting along the 
direction of the pulse propagation [1^ . 

We have analyzed the system m numerically and an¬ 
alytically; in Fig. 1, we display the intensity of the elec¬ 
tric field E in the left panels and the perturbation S' of 
the background radiation in the right panels. As an ini¬ 
tial condition, we use a Gaussian pulse for the electric 
field envelope, E = Aoexp {—+ {z + zo)^]/ao}, while 
the radiation perturbation is initially set to zero, and 
Vg zz c together with (3^ = 10“’^ (relevant to astrophysi¬ 
cal applications); see the upper panels of Fig. 1. The z 
derivatives in Eq. o are calculated numerically with a 
pseudo-spectral method, the r derivatives with a second- 
order difference scheme, and the system is advanced in 
time with a fourth-order Runge-Kutta scheme. 

Displayed in the middle panels of Fig. 1 is the initial 
pulse which has been moving in the z direction with a 
speed close to the speed of light. The pulse has self- 
focused, and exhibits a structure stretched along the z 
direction. The self-focusing can be understood in the 
framework of the two-dimensional nonlinear Schrodinger 
equation; during the initial phase of the pulse propaga¬ 
tion, the dispersive term plays a minor dynamical role. 
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FIG. 1: The collapse of the pulse and radiation gas distur¬ 
bance displayed using normalized electric field intensity \E'^ 
(left panels) and energy densities S (right panels). We have 
scaled the pulse electric field E by 2(3-^/reo/iC:)”^ and the 
radiation background S by 2Scl?>, while the time scale is nor¬ 
malized by (cfco)”^ and the length scale by . The initially 
Gaussian pulse has the amplitude Eq = ^eo/Sc)~^ 

and width uq = 50//co- The upper panels show the initial 
configuration, while the middle panels give the system af¬ 
ter 18 time units and the lower panels after 36 time units. 
The initial photonic pulse in the upper left panel self-focuses 
and then splits into two pulses, as can be seen in the left 
middle and lower panels. The energy background emanates 
Gherenkov-like radiation, clearly seen in the lower right panel 
as two yellow bands created behind the pulse. 


If the t and z variations are neglected in Eq. CB , we 
obtain S ~ 3//eo|-£^p. Thus, Eq. CB becomes the stan¬ 
dard, nonlinear Schrodinger equation in which the cubic 
nonlinear potential has the same sign as the diffraction 
term, supporting two-dimensional collapse of the pulse 
[ 3 . The time for complete two-dimensional collapse of a 
pulse can be estimated by means of Rayleigh-Ritz op¬ 
timization. Using normalized units (see Fig. 1), two- 
dimensional collapse will occur when the intensity and 
pulse width satisfies the inequality |Aop > l.Sttg The 
resulting collapse time tc — 4.3ao(AQao — 1.3)“^/^. In 
Fig. I, pulse splitting occurs after ^ tc/3, thus well be¬ 
fore field strengths reaches critical levels. 

However, here the collapse of the pulse is arrested due 


























3 


to the backscattering from the background; the core of 
the pulse slows down because of the nonlinear interac¬ 
tion with the background, while the flanks of the pulse 
continue to propagate with the speed of light, creating 
a fan-like structure in front of the slower moving pulse; 
see the middle left panel of Fig. 1. The pulse then splits 
into several parts with local maxima; one wider pulse 
which can be seen at k^z = 35 in the lower left panel fol¬ 
lowed by two smaller, narrower pulses at k^z « 27. We 
have carried out simulations with different values of /3^, 
which reveal that the flattening of the leading edge pulse 
is less pronounced for smaller values of j3z- The wider 
pulse is correlated with a slight depletion of the radia¬ 
tion energy density. Letting the pulse propagate further, 
the transverse variation for the wide pulse will become 
much smaller than the longitudinal variation, leading to 
an almost one-dimensional structure (depending weakly 
on the transverse coordinate r), thus making the nonzero 
dispersive term essential. This pulse moves with a super¬ 
sonic speed (> c/-\/3). For these broad supersonic pulses, 
we may solve Eqs. o for the stationary state. The acous¬ 
tic equation llbll gives ^ = —3/r(Ug-|-c^)/(3Ug —c^)eo|£'p. 
We note that the pulse will experience resonance phe¬ 
nomena as the group velocity approaches the sound 
speed. With the solution to Eq. 110, one can de¬ 
scribe the localized stationary solitary solutions for a 
pulse moving close to the speed of light according to 0 
\E\ ~ Eosech[(3£/5)^/^fco(2: —Ugt)], where Eq is the con¬ 
stant amplitude of the pulse, and e = cqEqIS’q is the 
relative energy density of the pulse. We point out that 
these pulses are not necessarily of extremely high ampli¬ 
tudes. The remnants of radially collapsing pulse asymp¬ 
totically form a train of low-amplitude axially modulated 
solitary pulses. Moreover, the pulses show a parabolic 
self-compression, which will compensate for the small, 
but non-zero, diffraction of the pulse, making the one¬ 
dimensional approximation valid over a longer period of 
time. In the final stage, we can see the formation of a 
‘photon wedge’ (seen in the lower right panel,) which is 
due to the pulse propagation with a supersonic speed, 
so that part of the energy of the pulse is released into 
Cherenkov-like radiation behind the pulse, analogous to 
the sonic shocks created behind supersonic flights in air. 
Thus, the dynamics of the photonic pulse is surprisingly 
complex, exhibiting a multitude of nonlinear phenomena. 

Simulations with higher initial amplitudes of the pulse 
show that the pulse again experience two-dimensional 
collapse in which the intensity may grow above the 
Schwinger limit before photon pulse splitting occurs. Be¬ 
cause of the nonlinear dominance, the speed of the pulse 
decreases below the sound speed of the background, and 
the Cherenkov energy loss through the photon wedges 
will therefore be small, reinforcing the photon pulse col¬ 
lapse. This mechanism leads to amplitudes where our 
model breaks down; thus, the photon collapse becomes 
catastrophic close to the Schwinger limit [l^. However, 


before the Schwinger limit is reached, higher order non¬ 
linear effects are likely to arrest the collapse [ 2 ^. This 
effect is similar to ultra-short intense laser pulses in air, 
where the formation of a plasma due to self-focusing gives 
rise to filamentation and halted collapse [23 . The plasma 
formation gives rise to a significantly longer propagation 
range for laser pulses in air, a behavior which the QED 
pulse propagation presented here is expected to share, 
but with the plasma formation replaced by higher order 
QED effects and the pair creation. 

Astrophysical environments can be of very extreme na¬ 
ture, exhibiting the largest energy levels known to man. 
In the case of regular neutron stars the surface magnetic 
field strengths reach 10^° — 10^^ G, while in magnetars 
they can reach 10^^ — 10^® G, the latter being close to the 
Schwinger limit. An interesting possibility arises in the 
context of neutron star and magnetar quakes, in which 
magnetic fields build up tensions in the star crust over 
long periods of time, and sudden bursts of energy are 
released from the star during the quakes. There, it is 
expected that large quantities of low-frequency photons 
would be ejected, forming an almost incoherent spec¬ 
trum of waves |29|. This photon gas could reach en¬ 
ergy densities tfo ~ 10^^ — 10^® Jm“^, corresponding to 
/i ~ 10“^° — 10“^. A short high intensity electromag¬ 
netic pulse, with wavelengths from the UV to gamma 
range (corresponding to 5 ~ 10“^° — 10“^), with its 
evolution modeled by the nonlinear Schrodinger equa¬ 
tion (Hall , passing through this low-frequency photon gas 
dynamically governed by the acoustic-like wave equation 
m, could be a source for gamma-ray bursts. The lat¬ 
ter are short (< tens of second) emissions of photons in 
the gamma range |^ . Apart from the pulse collapse 
described by the system o, there will also be a sig¬ 
nificant blue-shift of the photon pulse, due to the for¬ 
mation of steep intensity gradients in the system [T4|| . 
which we can see in the numerical results of Fig. 1. We 
note that the timescales for these events could be ex¬ 
tremely short depending on the frequency of the photons. 
One can therefore expect the formation of intense photon 
pulses with frequencies up to the gamma regime, within 
the incoherent photon gas created by magnetar quakes, 
possibly giving insight into the dynamics of gamma-ray 
bursts, since we here have a mechanism both for blue 
shifting, pulse compression and high intensity field gen¬ 
eration. However, we note that, in accordance with the 
standard relativistic fireball model, it would be more ap¬ 
propriate to model incoherent photons by a wave kinetic 
equation, instead of a Schrodinger equation. 

To summarize, we have considered the implications of 
the dispersive properties of the quantum vacuum for the 
case of intense photon pulses propagating on a radiation 
background. In the slowly varying acoustic wave limit, 
the pulse evolves similar to an ultra-short high-intensity 
pulse in nonlinear, normal dispersive media, with pulse 
collapse and splitting as a result. The analysis of the 
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full system of equations shows that the slowly varying 
acoustic limit is far from generic, and that the response 
of the radiation gas can have the same timescale as the 
pulse evolution. It is due to the self-generation of poten¬ 
tial wells, giving an attractive force between the photonic 
pulse peak and the acoustic disturbance. This can give 
rise to three-dimensional catastrophic photonic pulse col¬ 
lapse, where the pulse and radiation gas power increases 
towards the Schwinger limit. Moreover, given suitable 
initial conditions, the photonic pulses can evolve into a 
stable localized structure with high field strengths. The 
application of our work to astrophysical settings has been 
discussed. Specifically, the present investigation sheds 
light on gamma-ray burst dynamics, and gives a means 
for obtaining pulse intensities surpassing the ones achiev¬ 
able by known mechanisms. 
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